Physical Review E

Within the studies of stochastic gene expression, apart from the variability of copy number of gene products, the problems of threshold crossing of those products are biologically important as they often lead to terminal cellular events. Here, we study the threshold crossing problem of the messenger ribonucleic acid (mRNA) and present an exact probability distribution of first passage times in Laplace space. The function furnishes moments of any order and also predicts the characteristic time of the exponential tail of the distribution, which we match against Gillespie simulations. We find that all the measures of relative fluctuations of the threshold crossing times show U-shapes within this simple model of mRNA, as was found earlier in more mathematically involved models of threshold crossing time statistics of proteins. Furthermore, we extend the exact formula to include the phenomenon of DNA duplication and the corresponding doubling of transcription rate. As expected, the distribution varies considerably depending on the onset of the duplication stage within the cell cycle.

A fundamental biological process of transcription occurs in the cell nucleus, which is a complex medium that also contains multiple heterogeneous structures known as biomolecular condensates. Interestingly, some of these condensates contain mRNA molecules in addition to proteins, suggesting an important cellular role in transcription that is not yet well understood. In this work, we develop a minimal theoretical framework for quantitative investigation of the role of reversible mRNA condensation in transcription. Our discrete-state stochastic approach accounts for the most relevant processes, allowing us to explicitly evaluate the properties of the system and clarify the effects of condensation. Analytical calculations supported by computer simulations suggest that reversible mRNA condensation influences the transcription processes by maintaining a constant level of free mRNA in the nucleoplasm while lowering the degree of stochastic noise and increasing the robustness against external perturbations. Physicochemical arguments are presented to explain these observations. The proposed theoretical framework elucidates important microscopic aspects of transcription, providing a convenient quantitative tool for investigating complex biological phenomena.

Microbial populations generally evolve in fluctuating environments under time-varying conditions. These are often described by binary switching models, sometimes seen as coarse-grained feast-famine cycles, in which resource availability switches abruptly between abundant and scarce conditions. However, experimental studies suggest that feast-famine environments actually exhibit more complex temporal dynamics. Here, we study how two strains, one growing slightly slower than the other, compete for the same resources in fluctuating environments comprising a finite number of intermediate states, each having its own carrying capacity. Environmental switching between these states and their carrying capacities represents gradual changes in nutrient availability. This class of multi-state stochastic switching models can be interpreted as a coarse-grained description of feast-famine cycles and allows us to investigate strain competition under the gradual recovery and depletion of resources. By computational and analytical means, we characterise the population dynamics in these multi-state fluctuating environments. In particular, we study how the switching rates and distribution of carrying capacities affect the population-size statistics, fixation probability, and mean fixation time. By comparing these results with their counterparts in binary environments, we clarify how the frequency and amplitude of environmental fluctuations influence population dynamics in coarse-grained feast-famine cycles.

Efficient operation of neural networks has been linked to criticality in their underlying non-equilibrium excitation dynamics. However, obtaining experimental evidence of this conjecture remains challenging due to limited control and undersampling in biological systems. Here, we experimentally explore neural network criticality using an ultracold Rydberg gas as a highly controllable simulator. We highlight the similarity of the excitation spreading via Rydberg facilitation and the synaptic connection of spiking activity of neurons, giving rise to distinct absorbing and active phases. We systematically explore and resolve criticality criteria, including power-law scaling of excitation avalanches and the emergence of universal avalanche shape collapse. Crucially, we implement a controlled gain mechanism to compensate for atom loss, mimicking metabolic resource replenishment and stabilizing the system in a controlled non-equilibrium steady state. We find peak temporal correlations at the critical point and stochastic oscillations with dragon king avalanches in the active phase, consistent with predictions for systems orbiting criticality. Our work establishes facilitated Rydberg gases as a platform for investigating criticality, resource dynamics, and emergent oscillations in neural networks.

Based on our experimental observation of activation state oscillations of different frequencies used to communicate information by the master human stress response regulator protein p38 MAPK 1, we develop a simple graphical approach for understanding and predicting the behavior of complex biological networks acting upon signals carrying information as different frequency waves of chemicals. This approach uses the same techniques used for analyzing and understanding information transmission using waves of electrical currents and fields used in electrical alternating current (AC) circuits. We show how biological components can be organized to behave as standard components found in electronic telecommunications circuits. Finally, we demonstrate how such components can be organized into complex biological signaling cascades whose behavior can be qualitatively and quantitatively understood, and whose output resembles that experimentally observed in p38.

The study of tissue dynamics has been stimulated during the last decades thanks to the use of quantitative descriptions, with the development of several theoretical and computational frameworks, many of them revolving around the notion of reaction-diffusion systems, eventually with additional structure variables beyond time and space. The use of structure variables can accommodate phenotypic traits. In this work, we study a family of competition models, where a given population depends on a resource (e.g. oxygen) and several populations are competing for it. Our quantitative description incorporates phenotypic traits and heterogeneity at the level of cell cycle variations, which influence replication rates via oxygen consumption. This enables us to replicate the fitness of specific subpopulations to environmental conditions (e.g. oxygen shortage or external influences). Using numerical simulations, we show that such models display dynamical pattern formation in the form of coupled travelling wave profiles that expand or retreat at the same wave speed. The full theoretical analysis of such dynamics is quite involved; to circumvent this difficulty, we introduce a quasi-stationary approximation for the resource dynamics. We find that this approximation can reproduce the overall behaviour very accurately, with the additional benefit of allowing theoretical treatment of the reduced model. In this way, we provide estimates on the wave speed which are numerically shown to be robust across a wide range of macroscopic parameters of the full model. The wave speeds are thus found to depend strongly on the proliferation rate of the fittest population, resembling a winner-takes-all dynamics.

With 302 neurons and a rigorously characterized connectome, the nematode Caenorhabditis elegans represents a powerful model organism to study the fundamental roles of neuronal circuits in behavior. However, despite the breadth of research, many questions remain unanswered regarding how these organisms are able to successfully navigate their environment. Here, we present a biologically grounded dynamical circuit model for the investigation of sensory-guided behavior during C. elegans chemotaxis. Our mathematical model consists of the chemosensory neuron AWA, interneurons RIM and RIA, motor neurons, including SMDs and RMDs, and body wall muscles that provide proprioceptive feedback through stretch receptors. After optimization with an evolutionary algorithm, the model locomotes effectively toward a chemical attractant, realistically capturing nematode chemotactic behavior. Chemotaxis is ensured by sharp turns, which resemble the omega turns of living nematodes, as a key emergent property of the model. The sharp turning behavior is triggered by decreases in the concentration of the attractant. These result in reduced AWA activity, which in turn triggers disinhibition of RIM and subsequent changes in RIA oscillations. The ensuing coordinated changes in downstream motor neurons activity patterns produce sharp turns, which correct the nematodes path, so that the model worm heads toward the attractant, and remains at its proximity, after it reaches the gradient peak. The proposed framework, along with its emergent dynamics, provides new insights into the minimum requirements for C. elegans circuitry to display major features of its chemotactic behavior, including omega turns. In parallel, it generates experimentally testable hypotheses with respect to the participating neuronal elements.

Sensory integration within antennal lobe (AL) is thought to play a role in the guidance of odor tracking, though while responses to olfactory input within the AL have been well-studied, responses to mechanosensory stimuli, in the form of wind speed, have received less scrutiny. Recent experimental work has systematically characterized mechanosensory responses of individual neurons within the AL, showing four distinct response patterns, labeled as sustained, transient, biphasic, and offset types; furthermore, these experiments have demonstrated that the distribution of response patterns, as well as the response type of a fixed neuron, can vary with stimulus intensity (wind speed). In this work, we develop a realistic biophysical model of the AL and, using this model, we show that internal AL dynamics are capable of generating the response patterns observed experimentally - namely, we find that response type is determined by the interplay of slow synaptic inhibition, an intrinsic calcium dependent potassium (SK) current, and stimulus strength, and hence that a heterogeneous distribution of slow inhibition and SK current strength across the AL can lead to the emergence of all four types at a fixed wind speed. Moreover, similar to experiment, we find that the distribution of response patterns across the AL changes with wind speed. Finally, we examine the distribution of response types in the presence of a simulated odor stimulus as well as odor separation by the AL in the presence versus absence of a strong mechanosensory signal.

We show that the p53-Mdm2 gene interaction can be understood as a form of antithetic feedback, but with mixed feedback loops. In particular, we show that this viewpoint sheds light on the emergence and robustness properties of p53's oscillatory behaviour, as observed experimentally.

Maintaining a prescribed composition in engineered microbial consortia is difficult because small fitness differences can drive competitive exclusion. We study a two-strain consortium in continuous culture and develop a feedback architecture that regulates composition by selectively slowing the fast strain as a function of the population ratio. At the population level, we derive an idealized ratio-feedback law with a tunable positive coexistence equilibrium. We then propose a biomolecular realization using orthogonal quorum sensing, an sRNA-based ratiometric controller, and a ppGpp-mediated growth actuator. Exploiting the separation between slow population growth and faster intracellular controller dynamics, we use singular perturbation theory to show that, for sufficiently fast controller dynamics, the full implementation model inherits the coexistence equilibrium and its local stability properties from the reduced model. Numerical simulations validate the reduction and show how weaker timescale separation or loss of the assumed molecular regime degrades performance.

We investigate the dynamics of the spread of the infection due to a phage- during the range expansion of an E coli colony. We find a rich panoply of spatio-temporal dynamics at large colony-level length scales, arising due to the interplay of mechanical forces, growth and viral infection dynamics at smaller cellular length scales. Strikingly, the phages "surf" the front of the growing bacterial colony by hitchhiking on E coli that are advected due to the growth of bacteria in the colony, resulting in an anisotropic spread of the phages. We identify microscopic processes - of the phage release during lysis (occurring on short millisecond timescales at length scales close to that of a single bacterium) and local nematic alignment of the rod-like E coli bacteria (occurring on timescales comparable to the bacterial growth rate and the length scales of a few bacterial cells) - that enhance the advective effects driving the hitchhiking behavior. Combining our experiments with a mathematical model, we explain the multiple spatio-temporal dynamical regimes - from the coexistence of uninfected bacteria, resistant (lysogenic) cells and phages, to fixed points where the entire population turns resistant or remains uninfected. Our work broaches a new frontier in the investigation of eco-evolutionary processes arising from the mechanical and ecological forces due to the interactions between viruses and their hosts.

Cell motility is powered by the bacterial flagellar motor, a rotary nanomachine whose activity is dynamically modulated by external stimuli. In Magnetospirillum gryphiswaldense, magnetic and chemical inputs are thought to converge at the level of motility control, yet the mechanisms underlying chemotactic regulation remain poorly understood. Here, we seek to establish a quantitative reference-state of a model strain of magnetotactic bacteria through bacterial flagellar motor dynamics. To achieve this goal, magnetotaxis was kept as a natural factor by keeping the earth's magnetic field as the only source of magnetism. To avoid aerotaxis bias, oxygen gradients were removed by implementing two different constant oxygen conditions: environmental oxygen exposure, and limited oxygen exposure. Together, the tethered-cell bacterial flagellar rotational assay and the free-swimming assay in the absence of external stimuli presented in this paper establish a framework for investigating motor and cellular swimming adaptation to magnetic, aerotactic, and chemical signals. Under reference-state conditions, the bacterial flagellar motor showed a tendency to exhibit log-normal distributions for the time spent in each motor state: runs in different directions (counter-clockwise and clockwise), and pause. A semi-Markov reference-state model was developed to provide a quantitative description of bacterial flagellar motor dynamics. The model revealed that cell magnetic polarity modulates the transition pathways leading to the paused state, whereas oxygen, in the absence of a gradient, primarily regulates residence in each motility state.

Behavioral variability reflects differences in internal and external parameters and provides a substrate for natural selection. However, uncovering its structure is challenging because behavior evolves nonlinearly across multiple timescales. Our goal is to reconstruct phenotypic variation from finite behavioral observations without prescribing the scales at which individuals should be compared. To achieve this, we introduce a multiscale dissimilarity that compares dynamics across all statistically resolvable resolutions by constructing reduced transfer operators at progressively finer scales and weighting their differences by finite-sampling uncertainty. This dissimilarity defines a geometry of dynamical phenotypes, which we reconstruct using diffusion maps. We show that the resulting phenotypic coordinates locally follow the parameter directions to which dynamics are most sensitive. We validate our approach in a stochastic model and apply it to bacterial and larval zebrafish behavior. In bacteria, we identify run speed as a major component of inter-individual variability, associated with the chemotaxis protein CheB. In zebrafish, we find that prior prey experience reshapes behavioral dynamics and phenotypic variability, with paramecia experience concentrating fish in a prey-capture-associated region of phenotypic space. Together, these results provide a framework for identifying internal and environmental variables that structure dynamical phenotypes across scales.

We have studied the emergence of slow relaxation oscillations in next generation neural mass models with spike frequency adaptation. Relaxation oscillations connect low firing state (Down state) to high firing state (Up state) via the slow adaptation. In the examined cases, the orbit relaxes towards the Up State via a sequence of collective damped oscillations (peaks of activity), thus revealing population bursting dynamics. The slower is the adaptation time scale the higher is the complexity (number of peaks) displayed by the relaxation oscillations. In particular, a chaos-induced spike-adding mechanism regulates the increase in the number of peaks. In analogy to what found in the Hidmarsh-Rose neuron model, two different types of chaotic behaviors have been identified: Population Spiking and Population Bursting Chaos. The increase of the adaptation strength leads to shorter (longer) Up (Down) state durations somehow mimicking the effect of charbachol in in vitro experiments, where spontaneous slow waves are observed. Indeed, the scenario depicted in [1], where an increase of the concentration of carbachol induces a transition from anesthesia-like to sleep-like dynamics is consistent with our results based on the variation of the adaptation strength. HighlightsO_LISpike Frequency Adaptation (SFA) promotes the emergence of Slow Relaxation Oscillations C_LIO_LISpike-adding mechanisms, controlled by SFA, lead to Relaxation Oscillations of increasing complexity C_LIO_LITwo types of chaotic behaviours: Population Spiking and Population Bursting Chaos C_LIO_LISFA regulates Up and Down States durations and their correlation C_LI

Classical neurogeometric models describe the primary visual cortex as a fibered structure in which retinal position and local orientation are coupled through the geometry of the roto-translation group. We extend this approach to the visuomotor cortex by modeling it as an assemblage of visual and motor cortical geometries. The model combines orientation-selective representations, analogous to those of the primary visual cortex, with movement-direction-selective representations, analogous to those of the primary motor cortex, in order to describe the mixed visual and motor selectivity observed in the visuomotor cortex. We introduce a coupled visuomotor structure in which visual orientation and motor direction coexist over a common spatial plane and interact through a relative-orientation constraint. Neural responses are modeled by orientation- and direction-dependent profile functions, and preference maps are obtained from vectorized population responses. Numerical simulations generate visual, motor, and mixed visuomotor response maps. A competition rule between visual and motor responses produces incidence ratios close to experimental observations in macaque visuomotor cortex. This framework provides a first neurogeometric approximation of visuomotor functional architecture and a mathematical setting for studying visually guided action.

The relatively long flagella of spermatozoa from insects, birds, and octopuses display double waves, characterized by two superimposed helical waves. The prevalance of these highly organized waveforms across diverse taxa and distinct flagellar architectures hints at shared underlying physics, motivating a model of the flagellum as an elastic filament immersed in a viscous fluid, actively driven by a single set of internal bending moment waves. Simulations of a clamped filament show that it can buckle under its own activity into whirling and flapping states. A multiple-scales analysis of the elastohydrodynamic equations reveals how nonlinear interactions between fast undulations generate an effective compression driving buckling, and connects wave-driven buckling to classical follower-force instabilities. Extending the model to a swimming spermatozoon, the same instability produces double waves. Parameter estimates across species show that most observed double waves lie within the regime where buckling is permitted, supporting self-buckling as a generic physical mechanism for double waves.

In biological neural circuits, the dynamics of neurons and synapses are tightly coupled. We study the consequences of this coupling and show that it enables a novel form of working memory. In recurrent neural network models with ongoing Hebbian plasticity, we find that following oscillatory stimulation, neurons continue to oscillate long after the input is removed. This creates a dynamic form of memory that has no explicit storage or retrieval phases and that requires no prior knowledge of the input. We trace the mechanism of these "persistent oscillations" to an interaction between neurons and synapses that creates complex outlier eigenvalues of the connectivity matrix. This is shown both in simulation and analytically. We leverage this mechanistic understanding to generate persistent oscillations with prespecified dynamics, creating a dynamic analog of a classical Hopfield network. Our work demonstrates that coupling neuronal and synaptic dynamics enables novel forms of computation.

Recent advances in large-scale electrophysiological recording technologies now allow simultaneous measurement of spike trains from ensembles of neurons. A central challenge in mathematical neuroscience is therefore to identify functional assemblies and their collective organisation directly from this data. The goal of this work is to devise a method to infer the collective organisation of the neurons from their spiking activity. We construct weighted functional networks from neural spike trains using the van Rossum distance to quantify pairwise similarity. The functional network captures similarities in neuron firing patterns and provides a representation of their collective organisation. We hypothesise that similar neuronal assemblies will appear as clustered communities in the network, and employ the Louvain algorithm to detect such assemblies. We validate our approach using synthetic spike train data and simulated data from a stochastic block model of Leaky Integrate and Fire neurons subjected to external Poisson drives, where the ground truth is known. Finally, we apply our approach to large-scale recordings from the Allen Institute Visual Coding: Neuropixels dataset. We find that our methodology works well as long as a sufficient amount of data is available and the temporal structure is strong relative to the noise.

The dynamics of gene regulatory networks are governed by intrinsic noise, stemming from the random nature of biochemical reactions, and by extrinsic noise, arising from fluctuations in cellular components and environmental conditions. Together, these sources can compromise the reliability of predictive computational models if not properly accounted for, and capturing both effects within a single framework remains a non-trivial task in computational biology. In this work, we propose an uncertainty quantification framework that addresses these two contributions jointly: intrinsic stochasticity is described through a partial integro-differential equation (PIDE) for the protein probability density function, whereas extrinsic noise is represented as parametric uncertainty in the kinetic parameters. The propagation of the uncertainty is carried out via an intrusive polynomial chaos expansion (PCE), in which the PCE coefficients are obtained from a stochastic Galerkin projection of the PIDE, yielding a coupled deterministic system that is solved with standard numerical methods. We illustrate the approach on a positive autoregulatory gene network with one and two uncertain kinetic parameters. The proposed approach accurately reproduces the mean, variance, and full protein probability density function, including the bimodal distributions, at a substantially lower computational cost.

This study presents improvements to the Hodgkin-Huxley (HH) models of ionic conductance and action potential generation. Sodium and potassium conductances are expressed by a single analytical formula describing the impulse response of a convolution of exponential distributions within a short-memory integration space. Treating transmembrane ion transit duration as a random variable, conductance profiles are interpreted as realizations of the probability density functions governing ionic movements. Applying the central limit theorem, the lognormal distribution emerges as the asymptotic profile of ionic conductances, constituting a fundamental primitive for such biosignals. A temporal state-transition paradigm describes the action potential waveform through four successive membrane potential transitions. Applied to electrophysiological recordings from lamprey reticulospinal neurons, this framework enables indirect estimation of key physiological quantities, including depolarization threshold, Nernst potentials, and net ion fluxes across the membrane. These advances open new perspectives for parameter estimation from experimental data and neuronal network simulation.